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We investigate the noise properties of pump currents through molecular wires and coupled quantum dots. 
As a model we employ a two level system that is connected to electron reservoirs and is non-adiabatically 
driven. Concerning the electron-electron interaction, we focus on two limits: non-interacting electrons and 
strong Coulomb repulsion. While the former case is treated within a Floquet scattering formalism, we derive 
for the latter case a master equation formalism for the computation of the current and the zero-frequency 
noise. For a pump operated close to internal resonances, the differences between the non-interacting and 
the strongly interacting limit turn out to be surprisingly small. 
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1 Introduction 

Recent experiments with coherently coupled quantum dots [1^] and molecular wires [5, 6] deal with the 
transport properties of small systems with a discrete level structure. These experimental achievements 
generated new theoretical interest in the transport properties of nanoscale systems [7, 8]. One particular 
field of interest is the interplay of transport and electronic excitations by an oscillating gate voltage, a 
microwave field, or an infrared laser, respectively. Such excitations bear intriguing phenomena like photon- 
assisted tunnelling [3,9-16] and the suppression of both the dc current [17, 18] and the zero-frequency 
noise [19,20]. 

A further intriguing phenomenon in this context is electron pumping induced by a cyclic change of 
the parameters in the absence of any external bias voltage [21-23]. For adiabatically slow driving, the 
transfered charge per cycle is determined by the area enclosed in parameter space during the cyclic evo- 
lution [24,25]. This implies that the resulting current is proportional to the driving frequency and, thus, 
suggests that non-adiabatic electron pumping is more effective. For practical applications, it is also desir- 
able that the pump cuiTent flows at a sufficiently low noise level. It has been found that adiabatic pumps 
can be practically noiseless [26]. This happens, however, on the expense of having only a small or even 
vanishing current [27]. Outside the adiabatic regime, when the driving frequency is close to the internal 
resonances of the conductor, the cuiTent assumes much larger values while its noise nevertheless is clearly 
sub-Poissonian [28]. Since this prediction of an optimal working point has been made for non-interacting 
electrons, the question on the influence of Coulomb repulsion arises. 

An intuitive description of the electron transport through mesoscopic systems is provided by the Lan- 
dauer scattering formula [29, 30] and its various generalisations. In this formalism, both the average cur- 
rent [31] and the transport noise characteristics [30, 32] can be expressed in terms of the quantum trans- 
mission probabilities of the respective scattering channels. If one heuristically postulates that the current 
obeys a scattering formula, one should worry whether this complies with the Pauli principle or if it has to be 
ensured by introducing "blocking factors" [31]. For static conductors the current being the experimentally 
relevant quantity, is independent of these blocking factors, which renders this question rather academic. 
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This is no longer the case when the scattering potential is time-dependent. Then a scattered electron can 
absorb or emit energy quanta from the driving field, which opens inelastic transport channels [33-35]. So 
blocking factors indeed can have a net effect on the current, and it has been suggested to test the demand 
for them experimentally with driven conductors [36, 37]. In order to avoid such conflicts, one should start 
from a many-particle description. In this spirit, within a Green function approach, a formal solution for the 
current through a time-dependent conductor has been presented, e.g., in Refs. [36] and [38] without taking 
advantage of the full Floquet theory for the wire. Nevertheless in some special cases like, e.g., for conduc- 
tors consisting of only a single level [39,40] or for the scattering by a piecewise constant potential [41], an 
explicit solution becomes feasible. A complete Floquet theory provides in addition to a current formula a 
prescription for the computation of the Green function [16, 42]. 

The spectral density of the current fluctuations has been derived for the low-frequency ac conductance 
[43,44] and the scattering by a slowly time-dependent potential [45]. For arbitrary driving frequencies, the 
noise has been characterised by its zero-frequency component [19]. A remarkable feature of the current 
noise in the presence of time-dependent fields is its dependence on the phase of the transmission amplitudes 
[19,45]. By clear contrast, both the current in the driven case [19] and the noise in the static case [30] 
depend solely on transmission probabilities. 

When electron-electron interactions beyond the mean-field level become relevant, the direct application 
of a Landauer-like theory is no longer possible and one has to resort to other methods like, e.g., a master 
equation description for the reduced density operator of the wire [46^9]. For time-dependent conductors, 
this enables a rather efficient treatment of transport problems after decomposing the wire density operator 
into a Floquet basis. Then it is possible to study relatively large driven conductors [16] and to include 
also electron-electron interactions [50, 51] and electron-phonon interactions [52]. For the computation of 
the current fluctuations, one can employ a generalised master equation that resolves the number of the 
transported electrons. This degree of freedom is traced out after introducing a counting variable [53]. For 
various static transport problems, this approach has been followed by several groups [54-60]. 

After introducing our model, we review in Sec. 2 the Floquet scattering theory for the computation 
of the current and the zero-frequency noise. In Sec. 3, we derive a master equation approach which is 
applicable also in the presence of electron-electron interactions. These formalisms are in Sec. 4 employed 
for investigating the influence of Coulomb repulsion on the noise in non-adiabatic electron pumps. 



1.1 Wire-lead model 

A frequently used model for a nanoscale conductors like molecular wires or coupled quantum dots is 
sketched in Fig. 1 . It is described by the time-dependent Hamiltonian 

H{t) = i?wiro(i) + ^^loads + -f^contacts, (1) 

where the different terms correspond to the central conductor ("wire"), electron reservoirs ("leads"), and 
the wire-lead couplings, respectively. We focus on the regime of coherent quantum transport where the 
main physics at work occurs on the wire itself. In doing so, we neglect other possible influences originating 
from driving-induced hot electrons in the leads and dissipation on the wire. Then, the wire Hamiltonian 
reads in a tight-binding approximation with N orbitals |n) 

Hwiieit) — ^ ^ ^nn' {t)c\gC^i + -Hinteraction- (2) 

For a molecular wire, this constitutes the so-called Hiickel description where each site corresponds to one 
atom. The fermion operators c„s, c\^g annihilate and create, respectively, an electron with spin s = t, J. in 
the orbital |n). The influence of an applied ac field or an oscillating gate voltage with frequency f2 = 2'k/T 
results in a periodic time-dependence of the wire Hamiltonian: Hnn' = Hnn' {t)- For the interaction 

Hamiltonian, we assume a capacitor model, so that 

-ffintEraction — -^wire (-^wire 1) (3) 
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Fig. 1 Level structure of a double quantum dot with N = 2 
orbitals. The terminating sites are coupled to leads with chem- 
ical potential hl and /iir = //l + eV , respectively. 



describes the number of electrons on the wire. Below we shall focus on two 



where .?V„„-e = E„s cj, ^ 

limits, namely the interaction-free case U = and strong interaction, J7 — > oo, which finally means that 
the Coulomb repulsion is so strong that only states with zero or one excess electron play a role. 
The leads are modelled by ideal electron gases. 



■ff leads — ^liA^qs^'Lqs + 



(4) 



q,s 



where cj^^ (c^^) creates an electron in the state \Lq) (|Rg)) in the left (right) lead. The wire-lead turmelling 
Hamiltonian 



-^contacts = ^1 {y^l^^l^qs^^U + ^R9«4gaCjVs) + ll-C- 



(5) 



q,s 



establishes the contact between the sites \N) and the respective lead. This tunnelling coupling is 
described by the spectral density 



r,(e) = 27r^|y,,p5(e-e,) 



(6) 



of lead I = L, R. In the following, we restrict ourselves to the so-called wide-band limit in which the 
spectral density is assumed to be energy-independent, Tt{e) T^. 

To fully specify the dynamics, we choose as an initial condition for the left/right lead a grand-canonical 
electron ensemble at temperature T and electro-chemical potential /Ul/r- Thus, the initial density matrix 
reads 



-( if leads - /iL JVl - /liR JVr )/ fee T 



(7) 



where = cj^^c^^g is the number of electrons in lead f and k^T denotes the Boltzmann constant 



times temperature. An applied voltage V maps to a cheinical potential difference Hr — Hl = eV with 

being the electron charge. Then, at initial time to, the only nontrivial expectation values of the wire 
operators read (cJ,^,g,C£gs) = fe{eq)Su'5gg'6ss' where /^{e) = (1 + exp[(e - ij,i)/kBT])~^ denotes the 
Fermi function. 



1 .2 Charge, current, and current fluctuations 

To avoid the expUcit appearance of commutators in the definition of correlation functions, we perform the 
derivation of the central transport quantities in the Heisenberg picture. As a starting point we choose the 
operator 

Qt{t) = eNe{t) - eNi{to) (8) 

which describes the charge accumulated in lead £ with respect to the initial state. Due to total charge 
conservation, Q( equals the net charge transmitted across the contact i; its time derivative defines the 
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corresponding current 

m = ^^Qdt)- (9) 

The current noise is described by the symmetrised correlation function 

Si{t,t') = ^{[AIi{t),AIe{t')]+) (10) 

of the current fluctuation operator AIe{t) = Ie{t) — {Ie{t)), where the anticommutator S]+ = AB + 
BA ensures hermiticity. At long times, St{t, t') = Sf {t + T, t' + T) shares the time-periodicity of the 
driving [42]. Therefore, it is possible to characterise the noise level by the zero-frequency component of 
Se{t, t — t) averaged over the driving period, 

Se = ^ [ dt [ dTSe{t,t-T). (11) 

-< Jo J-oo 

Moreover for two-terminal devices, Se is independent of the contact i, i.e., Sl = <Sr = S. 

The evaluation of the zero-frequency noise S directly from its definition (11) can be tedious due to the 
explicit appearance of both times, t and t — r. This inconvenience can be circumvented by employing the 
relation 

^((O'W) - {Qi{t)f) = 2 1 dTSe{t,t-T), (12) 

which foUows from the integral representation of Eqs. (8) and (9), Qe{t) = J*^ dt' Ie{t'), in the limit 
to — > — oo. By averaging Eq. (12) over the driving period and using S{t, t — t) = S{t — T,t), we obtain 

S = {^^{AQj{t)))^, (13) 

where AQi = Qi — (Qe) denotes the charge fluctuation operator and (. . .)t the time average. The fact that 
the time average can be evaluated from the limit S = \uiita^-oo{AQ\ {t))/{t — to) > allows to interpret 
the zero-frequency noise as the "charge diffusion coefficient". As a dimensionless measure for the relative 
noise strength, we employ the so-called Fano factor [61] 

F = 4rr (14) 

e|/| 

which can provide information about the nature of the transport mechanism [30, 62]. Here, I denotes the 
time-average of the current expectation value {hit))- Historically, the zero-frequency noise (11) contains 
a factor 2, i.e. S' = 2S, resulting from a different definition of the Fourier transform. Then, the Fano factor 
is defined as F = S'/2e\I\. 

1.3 Full counting statistics 

A more complete picture of the current fluctuations beyond second order correlations is provided by the 
fuU counting statistics. It is determined by the moment generating function 

cP{x,t) = {e'^'^'^h (15) 
and aUows the direct computation of the fcth moment of the charge in the left lead via the relation 



(16) 

Copyright line will be provided by the publisher 



adp header will be provided by the publisher 



5 



Subtracting from the moments the trivial contributions that depend on a shift of the initial values, one 
obtains the cumulants. They are defined and generated via the so-called cumulant generating function 
In (f){x, t) which replaces (f) in Eq. (16) [63], so that the fcth cumulant reads 



a(ix)'= 



ln0(x,i) 



x=o 



(17) 



In a continuum linoit for the leads, both the moments and the cumulants diverge as a function of time, 
and one focusses on the rates at which these quantities change in the long-time limit. This establishes 
between the first two cumulants and I{t) and S{t) the relations 



m= -ie|-C(x,i) 



S{t) 



■c{x,t) 



(18) 
(19) 



For driven systems, these quantities are time-dependent even in the asymptotic limit and, thus, we char- 
acterise the transport by the corresponding averages over one driving period. Then expressions (18) and 
(19) become identical to the previously defined time averages / and S, respectively. Herein we restrict 
ourselves to the computation of the first and the second cumulant, despite the fact that also higher-order 
cumulants can be measured [64, 65]. 



2 Floquet scattering theory 

We now derive from the model described in Section 1.1 in the absence of electron-electron interactions 
expressions for both the current through the wire and the associated noise by solving the corresponding 
Heisenberg equations of motions. Since for = 0, the both spin directions contribute independently to 
the current, we ignore the spin index which means that we consider the current per spin projection. We 
start from the equations of motion for the annihilation operators in lead i, 

Clq — —jf-lqClq — -Vgq Cng, (20) 

which are straightforwardly integrated to read 

Ciqit) = c,,(to)e-'^^<'(*-*'')/'' - ^Veq '\t e-'^'"^ l''c^,{t - r), (21) 

where nn denotes the molecular wire site attached to lead I, i.e., tt-l = 1 and tzr = N . Inserting (21) into 
the Heisenberg equations for the wire operators yields in the asymptotic limit —oo 

Cneit) =-^Yl Hn„n'{t) C^' (t) - C„, (t) + ^e{t), (22) 

c„(t) =- ^^ifnn'(i)cn'(i), n = 2,...,N-l. (23) 

n' 

For a energy-dependent spectral density Ti = T({e), the dissipative part of the Heisenberg equation (22) 
acquires a memory kernel which complicates not only its solution but also the derivation of a current 
formula. For details, we refer the reader to Ref. [16]. 
The influence of the operator- valued Gaussian noise 

^^(*) = E ^i'i e-'^^"**-*"^/'^ cM (24) 

Q 
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is fully specified by the expectation values {^i{t)) = and 

i^Unm) = Su'^ j dee--(*-*')/V£(e) , (25) 

which follow directly from the definition (24) and the initial conditions (7). It is convenient to define the 
Fourier representation of the noise operator, ^^(e) = / dtexp(iet//i)^£(t) whose correlation function 

{^\{e) (e')) = fe{e) 6{e - e') dee' (26) 



follows directly from Eq. (25). 



2 . 1 Retarded Green function 

The equations of motion (22) and (23) represent a set of linear inhomogeneous equations and, thus, can be 
solved with the help of a retarded Green function G{t, t') = —{i/h) U [t, t')6{t — t'), which obeys 

(i/i^ - Hit) + ^r) G{t, t') = 5{t - t'), (27) 

where T{t) = |l)rL(t)(l| + |A^}rR(t)(A^| and Ti-it) is the one-particle Hamiltonian corresponding to 
Eq. (2). At this stage, it is important to note that the propagator of the homogeneous equations obeys 
U {t, t') = U{t + T,t' + T). Accordingly, the Fourier representation of the retarded Green function 

G{t,e) = -- dT e'^^/^C/(i, t-T) = G{t + T,e) = V e-^'=^*G(*=) (e) (28) 

^ -^0 k=-oo 

is also T-periodic in the time argument, so that it can be represented as a Fourier series. Physically, 
the Fourier coefficients G^*^^ (e) describe the propagation of an electron with initial energy e under the 
absorption (emission) of photons for fc > (fc < 0). In the Umiting case of a time-independent 
situation, all sideband contributions with k ^ vanish and G(t, e) becomes time-independent and identical 
to G(°)(e). From the definition (27) of the Green function and its Fourier representation (28), it can be 
shown that the solution of the Heisenberg equations (22), (23) reads 



(29) 



where we have defined Gn.mit, e) = {n\G{t, e)\ni). 

Below, we need for the elimination of back- scattering terms the relation 

G\t,e')-G{t,e)= (ih^ - e' + e + iT^G\t,e')G{t,e). (30) 

A proof starts from the definition of the Green function, Eq. (27). By Fourier transformation with respect 
to t', we obtain the relation 

(ih^^+e-n{t) + ^T)G{t,e) = l, (31) 

which we multiply by G^(t, e) from the left. The difference between the resulting expression and its 
hermitian adjoint with e and e' interchanged is relation (30). 
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2.2 Current through the driven molecular wire 

Owing to charge conservation, the (net) current flowing from lead I into the molecular wire is deter- 
mined by the negative time derivative of the charge in lead I. Thus, the current operator reads In = 
ie[i/(t), iV^J/ft, where Ni = J^q ^\q'^tq denotes the corresponding electron number and — e the electron 
charge. From Eqs. (21) and (24) then follows 

■p /"OO 

Ht) = ^ dT{cUt)ci(t-r)+cI(t-r)ci(t)}-e{cl(t)CL(t)+CiWciW}. (32) 

This operator- valued expression for the time-dependent current is a convenient starting point for the eval- 
uation of expectation values like the dc and ac current and the current noise. 



2.2. 1 Time-average current 

To obtain the current {I-L{t)), we insert the solution (29) of the Heisenberg equation into the current oper- 
ator (32) and use the expectation values (26). The resulting expression 

= ^I^r^Im / defe{ey'^"'G*u{t,e)Gu{t-T,e)+2eT^lm J deMe)Gn{t,e) (33) 



still contains back-scattering terms Gn and, thus, is not of a "scattering form". Indeed, bringing (33) 
into a form that resembles the static current formula requires some tedious algebra. Such a derivation has 
been presented for the linear conductance of time-independent systems [66], for tunnelling barriers [67] 
and mesoscopic conductors [68] in the static case for finite voltage, and for a wire consisting of levels 
that couple equally strong to both leads [38]. For the periodically time-dependent case in the absence of 
electron-electron interactions, such an expression has been derived only recently [19, 42]. 

Inserting the matrix element (1| . . . |1) of equation (30) eliminates the back- scattering terms and yields 
for the time-dependent current the expression 



(34) 



where qi,{t) denotes the charge oscillating between the left lead and the wire. Obviously, since <7l(0 
is time-periodic and bounded, its time derivative cannot contribute to the average current. The time- 
dependent current is determined by the time-dependent transmission 



/>CXD 

TLR(t,e) =rLrRRe / dTe''^"'Glj,{t,e)G^N{t - T,e). 
Jo 



(35) 



The corresponding expression for T^hit, e) follows from the replacement (L, 1) <-> {R, N). We empha- 
sise that expression (34) obeys the form of the current formula for a static conductor within a scattering 
formalism. In particular, consistent with Refs. [31, 36], no "Pauli blocking factors" (1 — fg) appear in our 
derivation. 

The dc current obtained from (34) by time-averaging can be written in an even more compact form if 
we insert for the Green function the Fourier representation (28). This results in 



(36) 



where 



=^L^R|GS(.)|^ 



T«(e) =rKrL|G^i(e) 



(37) 
(38) 
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denote the transmission probabilities for electrons from the right to the left lead and vice versa, respec- 
tively, with initial energy e and final energy e + kh^l, i.e., the probability for an scattering event under the 
absorption (emission) of photons if > (fc < 0). 

For a static situation, all contributions with k vanish and Tj^{e) = Tp^Ce). Therefore, it is 
possible to write the current (36) as a product of a single transmission T{e), which is independent of the 
direction, and the difference of the Fermi functions, /R(e) — fhi^). We emphasise that in the driven case 
this no longer holds true. 

2.2.2 Noise power 

In order to derive a related expression for the time-averaged current-current correlation function (1 1), we 
insert the current operator (32) and the solution (29) of the Heisenberg equations of motion. Then, we 
again employ relation (30) and the shorthand notation ek = e + khfl, so that we finally obtain 



2 

/R(e)/R(efe) 



+ TrFl I Yl Ti.G%-'^ (e.) [g['P (e)] * - iG[-^'^ (e,) | '/l (e)/R(efe) } ^^^^ 

k' 

+ same terms with the replacement {L,l) {R, N). 

2.2.3 Floquet decomposition in the wide-band limit 

Solving the equations of motion (27) for the Green function is equivalent to computing a complete set of 
solutions for the equation 

ih^^im) = (WwireW - m)), (40) 

which is hnear and possesses time-dependent, T-periodic coefficients. Thus, it is possible to construct a 
complete set of solutions with the Floquet ansatz 

IV-aW) =cxp[(-ie„/fi-7„)i]|u„(t)), (41) 
l"a(i)) =X^exp(-ifcm)|ua,fe). (42) 

k 

The so-called Floquet States \ua{t)) obey the time-periodicity ofHwireC^) and have been decomposed into 
a Fourier series. In a Hilbert space that is extended by a periodic time coordinate, the so-called Sambe 
space [69], they obey the Floquet eigenvalue equation [70, 71] 

(Wwire(i) - iS - ih^) Kit)) = (e« - i/i7a)|w«(t)). (43) 

Due to the Brillouin zone structure of the Floquet spectrum [69, 70, 72], it is sufficient to compute all 
eigenvalues of the first Brillouin zone, —%Q./2 < ea < hn/2. Since the operator on the l.h.s. of Eq. (43) is 
non-Hermitian, the eigenvalues — ihja are generally complex- valued and the (right) eigenvectors are not 
mutually orthogonal. Thus, to determine the propagator, we need to solve also the adjoint Floquet equation 
yielding again the same eigenvalues but providing the adjoint eigenvectors It can be shown that 

the Floquet states \ua{t)) together with the adjoint states form at equal times a complete bi- 

orthogonal basis: {u'^{t)\u(3{t)) = 5a/3 and l^a(*))(^a (*)l = 1- A proof requires to account for the 
time-periodicity of the Floquet states since the eigenvalue equation (43) holds in a Hilbert space extended 
by a periodic time coordinate [70, 73]. 
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Using the Floquet equation (43), it is straightforward to show that the propagator can be written as 

U{t,t') = ^e-'(^"/'^-'^«)(*-*')|n„(t))(«J(t')|, (44) 

a 

where the sum runs over all Floquet states within one Brillouin zone. Consequently, the Fourier coefficients 
of the Green function read 

i df f°° 

GW(e) = --/ — e"="* / dTe''^/''U{t,t-T) (45) 
Jo T Jq 

^ l^a,fc'+fc)«,fc/| 

In general, the Floquet equation (43) has to be solved numerically. In the zero temperature limit consid- 
ered here, the Fermi functions in the expressions for the average current (36) and the zero-frequency noise 
(39) become step functions. Therefore, the integrands are rational functions and the remaining energy 
integrals can be performed analytically. 



3 Master equation approach 

In the presence of electron-electron interactions, an exact treatment of the electron transport within a scat- 
tering theory is no longer possible and a master equation formalism can be an appropriate tool for the 
computation of currents [13, 16,46,51,52]. Recently, master equations have been established for the com- 
putation of current noise of various static conductors as well [53-60]. In the following, we develop such 
an approach for the case of periodically time-dependent conductors. 



3 . 1 Perturbation theory and reduced density operator 

We start our derivation of a master equation formalism from the Liouville-von Neumann equation ihlZ^t) = 
[H{t), TZ{t)] for the total density operator TZ{t). By standard techniques we obtain the exact equation of 
motion 

d ~ i ~ 1 f°° 

— TZ{t) = —-[Hy,ij.e-lea.ds{t),7l{0)] — ^ J dr[ilwire-leads, [-ffwire-leadsli-T, t), (47) 

where the tilde denotes the interaction picture with respect to the lead and the wire Hamiltonian, X{t,t') = 

is the propagator without the coupling. Below we will employ Floquet theory 
in order to obtain explicit expressions for these operators. 

As already discussed above, the moment generating function = (exp(ix7VL)) contains the full 
information about the counting statistics. For its expUcit computation, we define in the Hilbert space of the 
wire the operator 

T{X, t) = trieads{e'^'^^7e(i)}, (48) 

whose Umit X ^ obviously is the reduced density operator of the wire, .?-'(0, t) = p. After tracing out 
the wire degrees of freedom, becomes the moment generating function (f>{x) = tr^ire It will prove 
convenient to decompose into a Taylor series, 

fc=l 

where the coefficients J^k = triads (-^l^) provide direct access to the moments {N^) = tr^ire J^k- 

Copyright line will be provided by the publisher 



10 



F. J. Kaiser and S. Kohler: Shot noise in non-adiabatically driven nanoscale conductors 



Our strategy is now to derive from the master equation (47) for the full density operator an equation 
of motion for the J^k- For that purpose, we transform the master equation for TZ back to the Schrodinger 
picture and multiply it from the left by the operator exp(ixAfL)- By tracing out the lead degrees of freedom 
and using the conmiutation relations [A^l, V^] = V and [A^l, V^] = —V\ we obtain 

^^(X, t) = {C+ (e'^ - l)J+ + (e-''^ - l)J-}T{x, <)• (50) 

In order to achieve this compact notation, we have defined the superoperators J± and the time-dependent 
Liouville operator 

C{t)X= -^-[H^ir,{t),X] 

T-1 poo f- 

+ ^J^ dry de[e''^( - cic\Xfi^{e) + d\Xcifj^{e) - X~c\cjj^{e) + ciX~c\hie)) 

+ e-'^"( - Xciclhie) + clXhhie) - c\ciXh{e) + CiXc\Me)) 
+ same terms with the replacement 1,L ^ N,R, 

(51) 

which also determines the time-evolution of the reduced density operator, p = C{t) p. The tilde denotes 
the interaction picture operator c = c(t, t — t) and the Fermi function of lead ^, while fe = l — fg. The 
current operators 

J+{t)X = ^J^ drj de{e''^clXci + e-''^clXci) h{e), (52) 

J_{t)X = ^J^ dr J de{e''^^ciXc\ + e-''^ciXc\)h{e), (53) 

describe the tunnelling of an electron from the left lead to the wire and the opposite process, respectively. 
Note that these superoperators still contain a non-trivial time-dependence stemming from the interaction- 
picture representation of the creation and annihilation operators of wire electrons. 

3.2 Computation of moments and cumulants 

For computation of the current (18) and the zero-frequency noise (19), we generahse the approach of 
Ref. [56] to the time-dependent case. Since we restrict the noise characterisation to the Fano factor, it is 
sufficient to compute the long-time behaviour of the first and the second moment of the electron number 
in the left lead. This information is fully contained in the time-derivative of the operator up to second 
order in x, for which we obtain by Taylor expansion of the equation of motion (50) the hierarchy 

p = Cit)p, (54) 
^1 = C{t) Ti + {J+ (t) + J- (t)) p, (55) 
= C{t) + 2(j+(t) + J-{t))Tx + (j+it) - J-{t))p. (56) 

The first equation determines the time-evolution of the reduced density operator, which in the long-time 
limit becomes the stationary solution pQ{t). Note that for a driven system, it generally is time-dependent. 
Replacing in Eq. (55) p by po and using the fact that tr^u-e CX = for any operator X, we obtain the 
stationary current 

I{t) = -etlvmc^l = -etrwire(J'+ -I- J-)po{t). (57) 
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The dc current follows simply by averaging over one driving period and one ends up with the current 
formula ofRef. [74]. 

The computation of T\ (t) is hindered by the fact that the inverse of a Liouvillian generally does not 
exist. For static systems this is obvious from the fact that the stationary solution fulfils Cpo = 0, which 
implies that C is singular. This unfortunately also comphcates the computation of the second cumulant 
and we proceed in the following way: We start from Eq. (12) which relates the zero frequency noise to 
the charge fluctuation in the leads and write the time derivative of the first and the second moment of the 
electron number in the left lead by the operators ^1,2- From the equations of motion (55) and (56), we then 
find S = trwiie{2( Jlf + J- — — {J+ — J-)p}, where we again used the relation trwire = 0. 
An important observation is now that the first part of this expression vanishes for oc pQ, which can 
easily be demonstrated by inserting the current expectation value (57). Since po trwire acts as a projector 
onto the stationary solution po, we can define the "perpendicular" part 

J^i± = J^i — Po tr„ire Ti, (58) 

which fulfils the relation trwire J^i± = and obeys the equation of motion 

j^ii. = m + {j+{t) + J-{t) - I{t))po{t), (59) 

We will see below that in contrast to Ti, the long-time limit of the traceless !Fi± can be computed directly 
from the equation of motion (59). Upon inserting Eq. (58) into the equation of motion (56), we finally 
obtain for the still time-dependent "charge diffusion coefficient" the expression 

S{t) = tr„i,e {2{J+ + J-)Tx^ - (J+ - J-)po}, (60) 

whose time-average finally provides the Fano factor F = S/el. 



3.3 Floquet decomposition 

The remaining task is now to compute the stationary solutions po{t) and Ti±_ (t) from the time-dependent 
equations of motion (54) and (55). Like for the computation of the dc current in our previous work [74], 
we solve this problem within a Floquet treatment of the isolated wire, which provides a convenient repre- 
sentation of the electron creation and annihilation operators. 



3.3.1 Fermionic Floquet operators 

In the driven wire Hamiltonian (2), the single-particle contribution commutes with the interaction term and, 
thus, these two Hamiltonians possess a complete set of common many-particle eigenstates. Here we start 
by diagonaUsing the first part of the Hamiltonian which describes the single-particle dynamics determined 
by the time-periodic matrix elements Hnn' (t)- According to the Floquet theorem, the corresponding (single 
particle) Schrodinger equation possesses a complete solution of the form 

|*«W>=e-'^"*/^|¥'aW). (61) 
with the so-called quasienergies and the T-periodic Floquet states 

|^„(t)) = ^e-''="*|(/.„,fc). (62) 

k 

The Floquet states and the quasienergies are obtained by solving the eigenvalue problem 

( ^)n|ffw (t) {n I - ifi^) |¥'a W> = £« (*)) , (63) 
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whose solution allows one to construct via Slater determinants many-particle Floquet states. In analogy to 
the quasimomenta in Bloch theory for spatially periodic potentials, the quasienergies e<j come in classes 

^a,k = + kMl, A; G Z, (64) 

of which all members represent the same physical solution of the Schrodinger equation. Thus we can 
restrict ourselves to states within one Brillouin zone Uke for example < < Ml. 

For the numerical computation of the operators po and ^ix, it is essential to have an explicit expression 
for the interaction picture representation of the wire operators. It can be obtained from the (fermionic) 
Floquet creation and annihilation operators [16] defined via the transformation 

Cas{t) = ^{LPa{t)\n)Cns- (65) 
n 

The inverse transformation 

Cns = '^{n\iPa{t))Cas{t) (66) 

a 

follows from the mutual orthogonality and the completeness of the Floquet states at equal times [70]. 
Note that the right-hand side of Eq. (66) becomes time independent after the sunmiation. The Floquet 
annihilation operator (65) has the interaction picture representation 

Cas{t,t') = U^o{t,t')Cas{t)Uo{t,t') (67) 

= e-'(^«+^^-")(*-*')/'^c„,(i'), (68) 

with the important feature that the time difference t — t' enters only via the exponential prefactor. This 
allows us to evaluate the r-integration of the master equation (54) after a Floquet decomposition. Relation 
(68) can easily be shown by computing the time derivative with respect to t which by use of the Floquet 
equation (63) becomes 

^c„«(i,i') = -^{e^ + UN^)cc,s{t,t'). (69) 

Together with the initial condition Ca{t', t') = Ca{t') follows relation (68). Note that the time evolution 
induced by Wwire(t) conserves the number of electrons on the wire. 



3.3.2 Master equation and current formula 

In order to make use of the Floquet ansatz, we decompose the master equation (54) and the current formula 
(57) into the Floquet basis derived in the last subsection. For that purpose we use the fact that we are finally 
interested in the current at asymptotically large times in the limit of large interaction U. The latter has the 
consequence that only wire states with at most one excess electron play a role, so that the stationary density 
operator po (t) can be decomposed into the 2A'' -|- 1 dimensional basis { |0) , cl^g (t) |0) }, where |0) denotes 
the wire state in the absence of excess electrons and s =1, J,. Moreover, it can be shown that at large times, 
the density operator becomes diagonal in the electron number N^^e, so that a proper ansatz reads 

Po(t) = |0)/)oo(i)(0| + J2 cL|0)Pa.,/3.'(i)(0|c^.'. (70) 

a,/3,s,s' 

Note that we keep terms with a ^ (3, which means that we work beyond a rotating-wave approximation. 
Indeed in a non-equilibrium situation, the off-diagonal density matrix elements pap will not vanish and 
neglecting them might lead to artefacts [16, 75]. 
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By inserting the decomposition (70) into the master equation (54), we obtain an equation of motion 
for the matrix elements Pas,i3s' ~ {OlcasPviireC^ps' l^)- We evaluate the trace over the lead states and com- 
pute the matrix element {0\cas{t) . . . c^^, (i)|0). Thereby we neglect the two-particle terms which are of 
the structure cl^gcj^jO){0\ci3sCas- Formally, these terms drop out in the limit of strong Coulomb repul- 
sion because they are accompanied by a rapidly oscillating phase factor cxp(— it/A^wireT/S^)- Then the 
T-integration results in a factor fh{ea,k + U) which vanishes in the limit of large U. Since the total 
Hamiltonian (1) is diagonal in the spin index s, we find that the density matrix elements Pas,i3s' are spin- 
independent as well so that after a transient stage 

PaT,/3T(*) = Pal,l3l{t) = Pa/sit) (71) 

and Pat, 01 = 0. Moreover, the stationary density operator (70) obeys the time periodicity of the driving 
field [16] and, thus, can be decomposed into the Fourier series 

Pa/3(t) = ^e-''=^Va/3,fc (72) 
k 

and poo{t) accordingly. 

After some algebra, we arrive at a set of coupled equations of motion for Pa^it) which in Fourier 
representation read 

i(ea - e/3 - kMl)pal3,k =^ ^ {'Pa,k' +k"\^) {'^\<^l3,k+k") POO,k' {fl.{^a,k'+k") + /L(e/3,fc+fe" )) 

k'.k" 

Y ^ {'Pa,k'+k"\i){M'Pa',k+k")Pa'l3,k' f{£a',k+k") 

a' ,k\k" 

^ ^ (95/3%fe'+fc"|l)(l|¥'/3,fc+fc")PQ/3',fc' /(e/3',fc'+fc") 

6\k',k" 

+ same terms with the replacement 1, L ^ A'", R. 

(73) 

In order to solve these equations, we have to eliminate poo.fe which is most conveniently done by inserting 
the Fourier representation of the normalisation condition 

trwire PO {t) = Poo {t) + 2^^ Paa{t) = l- (74) 
a 

In order to obtain for the current an expression that is consistent with the restriction to one excess 
electron, we compute the expectation values in the current formula (57) with the reduced density operator 
(70) and insert the Floquet representation (66) of the wire operators. Performing an average over the 
driving period, we obtain for the dc current the expression [74] 

—^^^ (^^{'P0,k'+k\i){M'Pa,k)Pa0,k'fL{ea,k) 

— y~^(ya,fc'+fc|l}(l|ya,fc)p00,fc'/L(ea,fc, 
k' 

Physically, the second contribution of the current formula (75) describes the tunnelling of an electron from 
the left lead to the wire and, thus, is proportional to poofh which denotes the probabiUty that a lead state is 
occupied while the wire is empty. The first terms corresponds to the reversed process namely the turmeUing 
on an electron from site 1 1) to the left lead. 
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The decomposition of the equation of motion (59) for the long-time limit of T\\^ and the subsequent 
computation of the S from Eq. (60) proceeds along the same Unes with the only difference that the current 
operators J± yield an inhomogeneity and that the r.h.s. of the trace condition (74) is 

trwire ^1_L = (-?^1±)00 + 2 ^(J^l_L)aa = 0. (76) 

a 

3.4 Spinless electrons 

A particular consequence of strong Coulomb repulsion is the mutual blocking of different spin channels. 
This motivates us to also compare to the case of spinless electrons which is physically realised by spin 
polarisation. For spinless electrons, we drop in the Hamiltonian (1) all spin indices. Physically, this hmit 
is realised by a sufficiently strong magnetic field that polarises all electrons contributing to the transport. 
By the same calculation as above, we then obtain for the current also the expression (75) but without the 
prefactor 2. The factor 2 is also no longer present in the normalisation condition (74) which now reads 

trwire Po (i ) = Poo (i ) + X] Z^"" (*) = 1 ' ^''^ 

and accordingly in the equation of motion for T\^. 



4 Noise in non-adiabatic electron pumps 

In a mesoscopic conductor, a spatial asymmetry together with an ac driving can induce a pump current, 
i.e. a dc current that flows even in the absence of any net bias [21-23]. More precisely, the central neces- 
sary condition for this effect is the absence of a particular symmetry, namely generaUsed parity defined as 
the invariance under spatial reflection in combination with a time shift by half a driving period [16]. For 
adiabaticaUy slow driving, the transferred charge is determined by the area enclosed in parameter space 
during one cycle of the periodic time-evolution [24, 25]. This imphes that the resulting dc current, apart 
from non-adiabatic corrections, is proportional to the driving frequency. For time-reversal symmetric con- 
ductors, the parameters do not enclose a finite area and, thus, one finds for small driving frequencies that 
the pump current obeys 7 oc fi^ [28]. 

In both the absence and the presence of time-reversal symmetry, the current increases with the driving 
frequency, which suggests that pumping is more effective beyond the adiabatic regime. Non-adiabatic 
electron pumping is particularly interesting because at internal resonances of the central system the pump 
current can assume rather large values [3, 12, 13J. In the absence of electron-electron interactions, the 
pump current, in addition, exhibits resonance peaks with a remarkably low current noise [28] and the 
question arises whether this favourable property persists once Coulomb repulsion becomes relevant. We 
consider the setup sketched in Fig. 1 and described by the Hamiltonian (1) and compare the results for 
non-interacting electrons with those for strong Coulomb repulsion. Thereby we focus on the parameter 
set for which recently a low Fano factor has been predicted [28]: It is characterised by a large internal 
bias, intermediately strong wire-lead coupling, and resonant driving. Moreover at the prime resonance, 
the driving amplitude has to fulfil the condition J\ (A/ Wl) = -^/S/SFl.r, where J\ denotes the first-order 
Bessel function of the first kind. 

Figure 2a shows the current obtained for non-interacting electrons and for strong Coulomb repulsion 
with and without considering the spin degree of freedom. The dc current exhibits characteristic peaks 
whenever the fcth-order resonance condition = [(ei — e-^f" -\- A^]^/^/fcfi is fulfiUed. Interestingly 
enough, the magnitude of the current is practically the same for all three cases: The interaction yields 
at most differences of the order 10%. We also find that for U = oo, the resonance peaks are slightly 
sharper, in particular if one takes the electron spin into account. This behaviour has also been found for 
photo-assisted transport in molecular wires, where it is even significantly more pronounced [74]. The zero- 
frequency noise shown in Fig. 2b possesses a more involved double-peak structure with a local minimum 
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Fig. 2 DC current (a), zero-frequency noise (b), and 
Fano factor (c) for the non-adiabatically driven elec- 
tron pump sketched in Fig. 1 as a function of the driv- 
ing frequency. The dot levels possess the energies 
ei.2 = ±2.5A and are coupled to the leads with the 
strength Fl.r = 0.3A. The driving amplitude is 
A = 3.7 A and the temperature ksT = 0.005A. The 
peak at n w A/h corresponds to the first-order reso- 
nance, while the peaks at lower frequencies are higher- 
order resonances. 



n [A/h] 

at the centre. This results in a rather pronounced minimum of the Fano factor (Fig. 2c) which assumes 
values as low as F w 0.25, which means sub-Poissonian current noise. 

Our observations lead us to the conclusion that interactions are not an essential obstacle for tuning the 
double-dot electron pump into a low-noise regime as suggested in Ref. [28]. The reason for this is that for 
the pump configuration skechted in Fig. 1, one energy level lies below the Fermi energy while the other 
lies well above. Consequently in equilibrium for a sufficiently small wire-lead coupling, the left site is 
occupied while the right site is empty, whatever the interaction strength. Thus, the double dot is populated 
with only one electron so that interactions become irrelevant. Unless the driving amplitude is huge, this 
occupation is altered only slightly. Consequently interactions do not modify the current significantly. We 
emphasise that for strong dot-lead coupling F and finite interaction U, these arguments no longer hold true. 

In the experiment of Ref. [2], a typical inter-dot coupling is A = 50/ieV. Then, an internal bias eo = 5 A 
corresponds to the resonance frequency — 5A/fi,Ri 2nx 60GHz. In both interacting and non-interacting 
electrons, a wire-lead couphng F = 0.1 A results in an optimised pump current of the order 200pA with a 
Fano factor F « 0.25. 



5 Conclusions 

Within the present work, we were particularly interested in the question whether Coulomb repulsion would 
deteriorate the desirable noise properties of the pump current found for non-interacting electrons. As a 
central result, we found that very strong interactions do not alter the picture, but merely lead to slight 
modifications of the current and the zero-frequency noise, which are of the order of 10%. Therefore the 
essential conclusion of Ref. [28] holds also true in the strongly interacting case: Low-noise operation of a 
non-adiabatic electron pump requires a large internal bias in combination with a strong inter-dot coupling 
and resonant driving. The consequence is that a properly tuned driving ampUtude provides a relatively 
large pump current with clearly sub-Poissonian noise. Furthermore, we found that the influence of the 
electron spin is also not significant. 
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For non-interacting electrons, the transport problem can be treated within a Floquet scattering formal- 
ism, which relates the dc current to the transmission probabilities of elastic and inelastic transmission 
channels. The latter account for the possibility of photon emission and absorption by a scattered electron. 
The zero-frequency noise, by contrast, is not only determined by probabilities, but also depends on the 

scattering phases. 

As soon as electron-electron interactions play a significant role, a scattering approach becomes of lim- 
ited use. In the case of weak wire-lead coupling, the transport problem can then be treated within a master 
equation approach that is based on second order perturbation theory in the coupling. While previous mas- 
ter equation studies of driven transport are restricted to the computation of the dc current, we developed a 
method to also compute the noise properties. Our method is based on a time-dependent generalised master 
equation for the full counting statistics. Concerning its solution, the time-dependence of the Liouville oper- 
ator represents a particular challenge, because it renders the asymptotic state of the system time-dependent 
as well. We coped with this difficulty by decomposing the master equation into a proper Floquet basis. 

Concerning the noise of non-adiabatic electron pumps, still a couple of intriguing questions remain. 
Thus far, we revealed that the two limiting cases of vanishing and strong interaction allow for low-noise 
pumping, but the transition between these limits may nevertheless bear surprises. Thus studying current 
noise for a finite interaction strength is highly desirable. Moreover, albeit for ac driven conductors exper- 
imentally beyond the present state of the art, cumulants of higher order and eventually the full counting 
statistics can be of interest. In that spirit, our study can only be a first step towards a more complete char- 
acterisation of the electron transport in the presence of Coulomb interaction and time-dependent fields. 
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